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Abstract 

We reconsider the problem of the formation of a large- 
scale magnetic field in the accretion disks around 
black holes. In contrast with previous work we take 
into account the nonuniform vertical structure of the 
disk. The high electrical conductivity of the outer 
layers of the disk prevents the outward diffusion of 
the magnetic field. This implies a stationary state 
with a strong magnetic field in the inner parts of the 
accretion disk close to the black hole, 
keywords: accretion, accretion disks — magnetohy- 
drodynamics — black hole physics 

1 Introduction 

Early work on disk accretion to a black hole ar- 
gued that a large-scale magnetic field of, for ex- 
ample, the interstellar medium would be dragged 
inward and greatly compressed by the accreting 
plasma (Bisnovatyi-Kogan & Ruzmaikin 1974, 1976; 
Lovelace 1976). Subsequently, analytic models of 
the field advection and diffusion in a turbulent disk 
suggested, that the large-scale field diffuses out- 
ward rapidly (Lubow, Papaloizou, & Pringle 1994; 
Lovelace, Romanova, & Newman 1994) and prevents 



a significant amplification of the external poloidal 
field by electrical current in the accretion disk. This 
has led to the suggestion that special conditions (non- 
axisymmctry) are required for the field to be advected 
inward (Spruit & Uzdensky 2005). 

We reconsider the question of the advec- 
tion/diffusion of a large-scale magnetic field in a tur- 
bulent plasma accretion disk, taking into account its 
nonuniform vertical structure. The high electrical 
conductivity of the surface layers of the disk, where 
the turbulence is suppressed by the radiation flux and 
the relatively high magnetic field, prevents outward 
diffusion of the magnetic field. This leads in general 
to a strong magnetic field in the inner parts of accre- 
tion disks around black holes. 

2 The fully turbulent model 

There are two limiting accretion disk models which 
have analytic solutions for a large-scale magnetic field 
structure. The first was constructed by Bisnovatyi- 
Kogan and Ruzmaikin (1976) for a stationary non- 
rotating accretion disk. A stationary state in this 
disk (with a constant mass flux onto a black hole) 
is maintained by the balance between magnetic and 
gravitational forces, and thermal balance (local) is 



maintained by Ohmic heating and radiative radiative 
conductivity for an optically thick conditions. The 
mass flux to the black hole in the accretion disk is de- 
termined by the finite conductivity of the disk matter 
and the diffusion of matter across the large-scale mag- 
netic field as sketched in Figure 1. The value of the 
large-scale magnetic field in stationary conditions is 
determined by the accretion disk mass, which in turn 
is determined by the magnetic diffusivity of the mat- 
ter. For a laminar disk with Coulomb conductivity 
(which is very large), the mass of the stationary disk 
is also very large making the disk's self-gravity impor- 
tant. Correspondingly, the magnetic field needed to 
support a mechanical equilibrium is also very large, 
reaching in the central parts of the disk <~ 10 11 G 
for a stellar mass black hole, for a temperature and 
density at infinity T m <~ fO 4 K and for ~ 1CP 24 
g/cm 3 . The stationary magnetic field increases with 

3/2 

the black hole mass as ~ M fc ^ , and with a mass flux 

~ M 3/2 , where M ~ PooT^ 1 " 1 (Bisnovatyi-Kogan 
and Ruzmaikin 1976). 




Figure 1: Sketch of the poloidal magnetic field 
threading an accretion disk (from Bisnovatyi-Kogan 
& Ruzmaikin 1976). The field strength increases with 
decreasing radius owing to flux freezing in the accret- 
ing disk matter. 

It is widely accepted that the laminar disk is un- 



stable to different hydrodynamic, magnctohydrody- 
namic, and plasma instabilities which implies that 
the disk is turbulent. In X-ray binary systems the 
assumption about turbulent accretion disk is neces- 
sary for construction of a realistic models (Shakura 
and Sunyaev 1973). Therefore the turbulent accre- 
tion disks had been constructed also for non-rotating 
models with a large-scale magnetic field. A turbulent 
magnetic diffusivity was considered by Parker (1971), 
and by Bisnovatyi-Kogan and Ruzmaikin (1976). In 
the last paper the turbulent diffusivity was scaled by 
the parameters of the turbulent motion, similar to the 
scaling of the shear a- viscosity in turbulent accretion 
disk in binaries (Shakura & Sunyaev, 1973) where 
the viscous stress tensor component t r< p = aP, with 
a < 1 a dimensionless constant and P the pressure in 
the disk midplane. In a more consistent representa- 
tion, the coefficient of turbulent kinematic viscosity 
v in the Navier-Stokes equation is taken in the form 
v = (2/3)av s h, where v s = \fP~Jp is the isothermal 
sound speed and p is the midplane density of the disk. 
Using this representation, the expression for the tur- 
bulent electrical conductivity at is 



Mnh^/P/p 

Here, a — aia 2 . The characteristic turbulence scale 
is i — a\h, where h is the half-thickness of the disk, 
the characteristic turbulent velocity is v t = a-i y/P~fp. 
The mass of the turbulent magnetized disk is orders 
of magnitude less that in the laminar disk case. 

The evolution of a large-scale magnetic field 
threading a turbulent Keplerian disk can be esti- 
mated easily. This field arises from two sources: ex- 
ternal electrical currents and currents in the accre- 
tion disk. Evidently, the field generated by the cur- 
rents in the disk can be much larger than that due 
to the external currents. The magnetic field may be- 
come dynamically important, influencing the accre- 
tion disk structure and leading to powerful jet forma- 
tion, only if it is strongly amplified during the radial 
inflow of the disk matter. This amplification is possi- 
ble only when the radial accretion speed of matter in 
the disk is larger than the outward diffusion speed of 
the poloidal magnetic field due to the turbulent dif- 
fusivity rjt = c 2 /(inert)- Estimates by Lubow et al. 
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(1994) shown that for a turbulent conductivity (1), 
the outward diffusion speed is larger than the accre- 
tion speed. Thus it appears that there is no large- 
scale magnetic field amplification during Keplerian 
disk accretion. This conclusion is discouraging be- 
cause the most plausible models of many phenomena 
observed in the the systems with black holes inside 
the galaxy, as well as in the extragalactic supermas- 
sive black holes are connected with a large values of a 
large-scale magnetic fields. Nevertheless, this result 
directly follows from the equations of the standard 
disk structure, with the turbulent electric conductiv- 

ity (1) 

Lubow et al. (1994) did numerical calculations 
for a simplified situation with constant relative disk 
thickness h/r, constant kinematic viscosity v, and 
turbulent conductivity cr t . It is easy to show that 
the same result follows analytically for the standard 
accretion disk structure which can be written as 



M = Airpv r rh , h = v s /Qk 



4irr 2 haP = M(j-j m ) 



-n K aPh 
2 



2aT 4 c 



(2) 
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where vk = tQk is the Keperian velocity (e.g., 
Bisnovatyi-Kogan & Lovelace 2001). For regions far 
from the inner disk boundary, the specific angular 
momentum j = tvk ^> jin- The characteristic time 
t V i SC of the accretion disk matter advection due to the 
shear viscosity is t visc — r/v r . From the first three 
relations in equation (2), we obtain 



tvisc 

V r 



avt 



(3) 



To estimate the time-scale of outward magnetic field 
diffusion, we use 



tdiff = 



r\ r B r 



(4) 



(Lubow et al. 1994), where B r and B z are the large 
scale field components evaluated at the top surface of 
the disk. Here the coefficient of the magnetic turbu- 
lent diffusivity i] is obtained from (1) 
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(5) 



For stationary conditions, the large-scale magnetic 
field in the accretion disk is determined by the equal- 
ity tvisc = tdiff which implies 



~B~z 



a v s 
a vk 



a h 



a r 



< 1 . 



(6) 



In contrast, the coronal poloidal field solutions typi- 
cally have B r /B z ~ 1 at the disk surface (Bisnovatyi- 
Kogan and Blinnikov 1972; Ustyugova et al. 1999), 
which implies that tdiff *C t V i SC - This inequality in- 
dicates that the magnetic field is not amplified during 
accretion for these physical conditions. 

3 Turbulent disk with radiative 
outer zones 

Near the surface of the disk, in the region of low op- 
tical depth, the turbulent motion is suppressed by 
the radiative flux, similar to the suppression of the 
convection over the photospheres of stars with outer 
convective zones. The presence of the outer radiative 
layer does not affect the estimate of the characteris- 
tic time tvisc of the matter advection in the accretion 
disk because it is determined by the main turbulent 
part of the disk. The time of the field diffusion, on the 
contrary, is significantly changed, because the elec- 
trical current is concentrated in the radiative highly 
conductive regions, which generate the main part of 
the magnetic field. The structure of the magnetic 
field with outer radiative layers is shown schemati- 
cally in Figure 2. 

Inside the turbulent disk the electrical current is 
negligibly small so that the magnetic field there is al- 
most fully vertical, with B r <C B z , according to (6). 
In the outer radiative layer, the field diffusion is very 
small, so that matter advection is leading to strong 
magnetic field amplification. The field amplification 
will last until the magnetic forces in the region over 
the photosphere become of the order of the gravi- 
tational ones, and start to participate in the equi- 
librium balance. In such conditions the MHD and 
plasma instabilities are developed, decreasing the ef- 
fective electrical conductivity. We suppose, that in 
the stationary state the magnetic forces could sup- 
port the optically thin regions against gravity. In the 
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Figure 2: Sketch of the large-scale poloidal magnetic 
field threading a rotating turbulent accretion disk 
with a radiative outer boundary layer. The toroidal 
current flows mainly in the highly conductive radia- 
tive layers. The large-scale (average) field in the tur- 
bulent region is almost vertical. 

nonrotating magnetized disk magnetic forces support 
the whole disk against the gravity, so they should be 
much higher. When the magnetic force balances the 
gravitational force on the outer optically thin part of 
the disk of surface density T, p h one finds the following 
relation takes place 



GMY^ ph B Z I^ 



2c 



4"7T 



(7) 



(Bisnovatyi-Kogan & Ruzmaikin 1976). The sur- 
face density over the photosphere corresponds to a 
layer with effective optical depth close to 2/3 (e.g., 
Bisnovatyi-Kogan 2001). Wc estimate the lower limit 
of the magnetic field strength, taking n es (instead of 
the effective opacity K e // = ^K es n a ). Writing 



->ph 



2/3 



(8) 



we obtain = 5/3 (g/cm 2 ) for the opacity of the 
Thomson scattering, K es — 0.4 cm 2 /g. The absorp- 
tion opacity n a is much less than n es in the inner 
regions of a luminous accretion disk. Thus using in 
equation (7) from equation (8), we estimate the 
lower bound on the large-scale magnetic field of a 
Kcplcrian accretion disk as 



B z 



1 



10 8 G- 



1 



xwm 



(9) 



where x — r/r g and 



M/Mq. For comparison, 



the surface density of the disk in the inner ra- 
diation dominated region, where we may expect the 



largest values of the magnetic fields, is 



_ 80V2x 3 / 2 / 
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(10) 



(see Bisnovatyi-Kogan 2001). The maximum mag- 
netic field is reached when the outward magnetic 
force balances the gravitational force on the disk 
of surface mass density E p ^. In equilibrium, B z ~ 
^/S p h- We find that B z in a Keplerian accretion 
disk is about 20 times less than its maximum possi- 
ble value for x =10, a = 0.1, and m = 10. 



4 Discussion 

An important question is the energy density of the 
large-scale stationary magnetic field in comparison 
with the rotational or gravitational energy density of 
the disk. For a non-rotating magnetized accretion 
disk, the energy density of the field is of the order 
of the gravitational one so that the magnetic field 
strength is very large in the vicinity of a black hole 
(Bisnovatyi-Kogan & Ruzmaikin 1976), and may be 
many order of magnitude larger than the external 
seed field. In the case of the fully turbulent Keplerian 
disk the poloidal magnetic field tends to drift outward 
(Lubow et al. 1994; Lovelace et al. 1994) so that its 
value cannot significantly exceed the strength of the 
large-scale seed magnetic field. 

The suggestion of a fully turbulent accretion disk 
with a small turbulent conductivity is violated in 
the outer surface layer of the accretion disk where 
the optical depth is small and the turbulence is sup- 
pressed by the strong radiative flux. This is simi- 
lar to the radiative layer above the main body of a 
convective star (similar to the Sun; e.g., Bisnovatyi- 
Kogan 2001). In the radiative layer, the magnetic 
field diffusion is much slower than in the region of 
fully developed turbulence. In the radiative layer the 
diffusion is determined by the classical Coulomb con- 
ductivity which is very large. The diffusion in this 
layer is practically negligible. The electric current 
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is concentrated in the radiative layer, and the main 
body of the turbulent disk is almost current-free and 
thus force- free. The magnetic field lines in this region 
are almost straight as shown in Figure 2. Because 
of the negligible diffusion in the radiative layer, the 
large-scale field drifts inward until the dynamical ac- 
tion of the magnetic field on the photosphere becomes 
comparable with that of centrifugal and gravitational 
forces. At this point the inward drift of the field will 
be halted and a the stationary state formed where the 
the magnetic, centrifugal, and gravitational forces on 
the optically thin region will be comparable, and de- 
viations from the Keplerian angular velocity will be 
significant. The strength of the magnetic field for 
such conditions is smaller than in the nonrotating 
disk of Bisnovatyi-Kogan and Ruzmaikin (1976), but 
still it is very large in the vicinity of the black hole. 
In this situation we may expect a nonuniform distri- 
bution of the angular velocity over the disk thickness: 
The main body of the turbulent disk is rotates with 
the velocity close to the Keplerian one, and outer op- 
tically thin layers rotate substantially slower. 



Self-consistent models of the rotating accretion 
disks with a large-scale magnetic field requires so- 
lution the equations of magnetohydrodynamics. We 
expect two different self-consistent solutions for the 
same set of the input parameters. In the case of a 
fully turbulent disk without radiative surface layers, 
the large-scale magnetic field will remain close to the 
value of the seed field and the disk's angular velocity 
will be close to Keplerian. In the second solution the 
strength of the magnetic field is large, and it may 
greatly exceed the strength of the seed field. In this 
solution the angular velocity distribution may deviate 
considerably from the Keplerian law. In presence of 
the radiative layers the solution with a small field will 
not be stationary, and a transition to the strong field 
solution will take place. We conclude therefore that 
the strong field solution is the only stable stationary 
solution for a rotating accretion disk. Further inves- 
tigation of the build up of strong large-scale fields by 
accretion is underway (Rothstein et al. 2007). 
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